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This is an updated review of recent work done by the authors on a proposal for non- 
perturbatively stable 2D quantum gravity coupled to c < 1 matter, based on the flows of 
the (generalised) KdV hierarchy [|T], ^. 

Since the discovery of the continuum limit of 2D quantum gravity coupled to minimal 
matter Q to all orders in the genus expansion, the question of a non-perturbative 
definition of these theories has remained. In refs.JI], |^ we have presented a 'physical' non- 
perturbative definition by elevating the fundamental (generalised) KdV flow symmetry, 
known to exist perturbatively, to a non-perturbative principle. Concentrating on KdV for 
simplicity, these flows are 

du , , , 

0^^= Rk+im = ik+iu (1) 

where 

Y.{k + l/2)tkRk -z = n = 0. (2) 

k 

Here the string susceptibility u = V" , prime is d/dz where z = ^/u, while /i and u are 
the renormalised 'cosmological constant' and The parameters t^/^ couple to the 

matrix-model integrated local operators Ok on the surfaces and Rk are the Gelfand-Dikii 
differentials The flows determine correlators of the local operators: 

92 d d 

__<0,...0.,>^5^...^. (3) 

The string equation (0) is the usual condition imposed by the hermitian matrix model 
(HMM) 1^, leading to non-perturbative instability in particular for pure gravity and 
probably all unitary c < 1 matter. Let us relax this condition and ask for the most 
general string equation compatible with the flows (|I|). To derive it note that dimensionally 
if [u] = 1 then [z] = —1/2 and [t^] = — (A; + 1/2) as follows from (0). Then u as a function 
of its dimensionful arguments must have the scaling symmetry^ 

J2ik + 1/2) tkj- + -zu + u = (4) 
k=o ^^fe ^ 

Using (|l]) and the recurrence relation [^]; 

R'k+i = \Rk - uR'u - \u'Rk (5) 

this becomes: 

^7^'" - uTZ' - ^u'n = (6) 
Multiplying by 71 and integrating gives the string equation [|^: 

un^ - ]^nn" + ^ (7^')' = o (7) 



^We assume initially that no new dimensionful parameter arises at the non-perturbative level. 
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The constant of integration (which is an order term) has been set to zero by hand 
in order that the HMM equation 7?. = satisfies (|^) asymptotically. This is necessary 
to reproduce the perturbative (genus) expansion at 2; — +00. Note that the asymptotic 
expansion of the solution to (0) is uniquely determined once we have fixed the spherical 
approximation. 

As an illustrative case we can consider pure gravity t^a = tm<^m2- In the spherical 
approximation, where we can neglect derivatives, (0 ) becomes: 

u{u^ -zf = (8) 

As z ^ +00 we take u = -y/i and the corrections will yield standard genus expansion, 
by construction. If we require u to have a real asjTiiptotic expansion at z ^ — 00 then 
we must take the root u = 0. The susceptibility then begins at the 'torus' level and in 
fact has non-vanishing contributions only every 2nd genus. More generally it is every m}^ 
genufl g: 

u = --^T.Cnl-^) ; Co = l/A ; .^-00 (9) 

^ n=0 ^ 

It is highly likely that all other non-perturbative solutions to (0) reproducing the standard 
genus expansion at +00 are complex or have poles on the real axis and are thus physically 
unacceptable[0]. The z ±00 behaviour is now enough to fix completely the boundary 
conditions for (^ and the (unique) numerical solution we found for pure gravity is 
shown in figure 1. On the question of uniqueness we note also that at tm = tm^mi, (0) 
is reducible to Painleve II in x by transformation u = 2%^ + z. With our physical 
boundary conditions it has already been shown analytically that there is a unique, real. 



pole-free solution |10 



Thus, starting from the KdV fiows and requiring physically acceptable behaviour in 
choosing boundary conditions, we are led to the solution for pure gravity shown in figure 
1. Considering u as the potential for a hamiltonian operator H, whose eigenvalues govern 
the non-perturbative positions of charges in a Dyson gas, we see that the spectrum is 
continuous and bounded below. In fact our analysis is equivalent to taking the continuum 
limit of an appropriate (critical) Dyson gas on the positive real line, as we show in a 
moment. We can compare our stabilisation with results provided by the SUSY D=l / 
stochastic quantisation of the HMM investigated in refs.p]. The non-perturbative differ- 
ence between the two at positive z is obvious from the manifestly different behaviour of 
the solutions at 2; — > — cxo. For the model investigated in [Q, since u is real and oc at 
z ^ —00 [|lT], |12|, the susceptibility cannot satisfy (|^) and so violates the KdV flows non- 
perturbatively. Also there still appear to be non-perturbative ambiguities in the approach 
of in its most general form, while there are none for our approach. 

Our analysis so far has dealt directly with the continuum limit and it is natural to ask 
whether there is an appropriate formulation in terms of a Dyson gas. The universality class 
of the critical behaviour in such a gas can be conveniently characterised by the structure 



*The same asymptotic behaviour has also recently been found in the unitary matrix model with 
'external fields' @ 
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of the eigenvalue density in the scaled neighbourhood of the end of an arc of eigenvalues, 



in the spherical approximation [|13| . Figure 2(a) shows this region as appropriate for pure 
gravity (m = 2) at z > 0. Eigenvalues are concentrated on the cut (wavy line) ended by 
a square root branch point, and formally extending the expression for the density off this 
cut, its first integral has the interpretation of effective potential for one eigenvalue, Ves 



T^. For the mth critical point there are m — 1 extra zeros of the density in the scaled 
neighbourhood of the branch point, shown as dots in the figure. The behaviour as a 
function of z is easily calculated by using the WKB approximation for the spectrum of the 
co-ordinate operator H, with the appropriate 'potential' u from (H). This also facilitates 
a comparison of the extremal values of V^s with the leading exponential corrections to the 
asymptotic solutions to (0) [0]. As 2; — > the branch point and the extra zero in figure 
2(a) collide and at first it is difficult to see how a real, sensible answer is maintained 
for z < 0. A solution is shown in figure 2(b) whereby the density now diverges, rather 
than vanishes, like a square root at the branch point, thus generating an extra zero. This 
scenario is thoroughly natural if there is a 'wall' in the problem, which the branch point 
hits at 2; = 0. In this case the two zeros depart into the complex plane perpendicular to 
the real axis as z becomes more negative. The structure shown in figure 2(a) to the left of 
this wall is then in the sense of analytic continuation (V^fj = 00 strictly). Simply speaking 
the wall will prevent eigenvalues leaking out of the arc and this is in fact the origin of 
non-perturbative stability in our formulation, as we now show. 

We can derive the full continuum limit a la Douglas which also indicates the inclu- 
sion of the generalised KdV hierarchy. Insertions of the eigenvalue operators —A, —d/dX 
have, for the mth critical point of the one-matrix model, the continuum limits: 

2m- 1 

-H = Q = d^-u ; P = ^id' ; ^ = ^/^^ (10) 

1=0 

The eigenvalue space of the HMM is IR and the canonical momentum operator is P 
(translations). The c.c.r. [P,Q] = 1 determines the a, : i < 2m — 1 uniquely (a2m-i 
is an undetermined non-universal constant), the order d'^ of this relation being TZ = 0. 
Our stabilisation scheme is equivalent to the restriction to 1R+ i.e. the imposition of a 
'wall' in the scaling region. To show this we note that the canonical momentum on 1R+ 
is P, the generator of scale transformations, satisfying the new c.c.r. [P, Q] = Q. HQ 
is conjugate to macroscopic loop length then these local scale transformations are related 
to physical scale transformations of the string/one-dimensional universe. We can find P 
by introducing the fractional-power pseudo-differential operator and its differential part: 

Qm+1/2 ^ ^2m+l ^^^^ ^/2){m, d^"^-'} + ■■■ + {Rm+l[u],d-'} ■■■ (11) 

[QT'^\Q] = li^^^, (12) 

We expect P = J2i=o~^ otid^ but Q™"*"^^^ does not quite do the job. To obtain a term d"^ 
on the r.h.s. of the c.c.r. we must use most generally: 

P = a,^+iQT"' - (13) 
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giving the (differential of the) string equation; 



[P, Q] = a2m+iC+i + — + £ = £ -u. (14) 

This is (^ at the mth critical point, modulo a2m+i which may be absorbed into the string 
coupling. By explicitly taking the continuum limit of the Dyson gas on 1R_|_, realised for 
example by the complex matrix model one finds that the constant arising from an 
integration of ([I^) is zero as in (|^). Moreover the physical boundary conditions we chose 
earlier are precisely the ones appropriate to this gas. In the spherical approximation u is 
the scaling part of the end of the eigenvalue density. For 2; > the 'wall' has no affect on 
the WKB expansion, in particular the leading term is y/z as in the HMM. At z < the 
end of the density remains fixed at the wall so u = in the spherical approximation. 

It is natural to suppose that these considerations have analogues for the generalised 
KdV hierarchy, related to the full set of (p, q) minimal models 0, H]. One would generalise 
to: 

q-2 

Q = d'^ + Y^ u,d' P = Q^'i (15) 
j=0 



Q^+^/^ - -^d (16) 



Some details of this will be explored in [|T5|. 



Let us finally consider the Dyson- Schwinger equations^ of the one-matrix models, at 
least that part of them which may be represented as formal Virasoro constraints [0. 
Our differentiated string equation (|^) follows from the constraint Lqt = (where u = 
—2(P logr and z ^ z + 1^) on using the KdV flow ([I|). This is not suprising in that it is 
precisely the expression of scale invariance. The higher Virasoro constraints: 

L„r = : n>l (17) 

where 

follow by applying the recursion operator which generates higher symmetries of the KdV 
hierarchy [00. The HMM string equation derives from L_it = and one might ask 
what happens to this constraint in our formulation. Since we are working on 1R_|_ there 
is no translation invariance and the L_i constraint is naively absent. More precisely the 
relevant DS equation picks up a boundary term from the 'wall': 

L-.r=|l (19) 

where a is the scaled position of the 'wall', which up to now we have taken as the origin. 
In fact there is no reason why we should choose a = and more correctly we might 

^Some of what follows was elaborated while these notes were being prepared. 

^H.Kawai informed me that the DS constraints {L^t = : n > 0} should in turn imply that r is a 
T- function of the KdV hierarchy subject to Li^t = (see his contribution to these proceedings). 
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consider it as an extra (non-perturbative) parameter u = u{z, tk, cr). The correct canonical 
momentum gets shifted now to P^- = P + aP and thus 

[Pa,Q + o-]=Q + a (20) 

implying 

^7^'" - uW - ^u'n + (xn' = (21) 

which may now be integrated to give the string equation. But cr, which has the same 
dimension as u, now contributes to the scaling equation (§) a term adu/da. Comparing 
with pl| ) we identify 

which is equivalent to (|T9|). The other Virasoro constraints also pick up boundary terms 
at general a 

Lnr = a-^'^ (23) 
da 

as follows simply from varying the boundary of the eigenvalue integration as cr ^ cr+ecr"+^. 
The parameter a can be set to zero by an analytic redefinition of the together with 
u — cx—^u. In this sense it is redundant and in fact if we assume, following the HMM, that 
the non-perturbative loop expectation is \I^(/) =< e^'"^ >, a corresponds to a boundary 
cosmological constant |[T^. The loop wavefunction in the scaling limit is 

^(/) oc / p{iP)e-^^dilj (24) 

where _ 

POO 

p{'4>) oc / < X I 5{'4> - H) \ X > dx (25) 
J fj, 

Hence for a > the wavefunction is guaranteed to have sensible exponentially decreasing 
large / behaviour. 
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